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Lattice-QCD-based Schwinger-Dyson Approach for Chiral Symmetry 
Restoration at Finite Temperature 

Hideaki lida'*, Makoto Oka'' and Hideo Suganuma'* 

''Dept. of Physics, Tokyo Institute of Technology, Ohokayama 2-12-1, Meguro, Tokyo 152-8551, Japan 

We propose the Schwinger-Dyson (SD) formaUsm based on lattice QCD, i.e., LQCD-based SD formalism, 
for the study of dynamical chiral-symmetry breaking in QCD. We extract the kernel function K{p^) in the SD 
equation from the lattice data of the quark propagator in the Landau gauge. As remarkable features, we find 
infrared vanishing and intermediate enhancement of the kernel function K{jP) in the SD equation. We apply the 
LQCD-based SD equation to thermal QCD, and calculate the quark mass function at the finite temperature. We 
find chiral symmetry restoration at the critical temperature Tc ~ lOOMeV. 



1. Introduction 

Dynamical chiral-symmetry breaking (DCSB) 
is one of the most important nonperturbative 
features in QCD. For the study of DCSB, the 
Schwinger-Dyson (SD) formalism |1I2I3I4| has 
been used as an interesting and powerful method. 
For QCD, the SD formalism consists of an infinite 
series of nonlinear integral equations which deter- 
mine the n-point Green's function of quarks and 
gluons, and therefore it includes the infinite-order 
effect of the QCD gauge coupling constant g. 

For instance, the SD equation for the quark 
propagator S{p) is described with the nonpertur- 
bative gluon propagator Dp,u{-p) and the nonper- 
turbative quark-gluon vertex gT,j{p, q) as 

S-\p) - 5o-i [ ^^S{q)D^,{p - q)T,{p, g), (1) 



where S'o(p) denotes the bare quark propagator 
and the simple notation = J been 
used in the Euclidean metric. 

In the practical calculation for QCD, however, 
the SD formalism is drastically truncated: the 
perturbative gluon propagator and the one-loop 
running coupling are used instead of the non- 
perturbative quantities in the original formalism. 
This simplification seems rather dangerous be- 
cause some of nonperturbative-QCD effects are 
neglected. 

In this paper, we formulate the SD equation 



based on the recent lattice QCD (LQCD) results, 

1. e., the LQCD-based SD equation, and aim to 
construct a useful and reliable analytic frame- 
work including the proper nonperturbative effect 
in QCD. Using the LQCD-based SD equation, we 
investigate also DCSB at finite temperatures. 

2. The Quark Propagator in Lattice QCD 

First, we briefly review the quark propagator 
S{p) in lattice QCD. The inverse quark propaga- 
tor in the Landau gauge is generally given by 



S-'{p)^Z{p')/{i> + M{p')} 



(2) 



in the Euclidean metric. Here, M{p'^) is called as 
the quark mass function, and Z{p'^) corresponds 
to the wave-function renormalization of the quark 
field. In the quark propagator, DCSB is charac- 
terized by the mass generation as ^ 0. 

The quark mass function MijP') in the Landau 
gauge is recently measured in lattice QCD at the 
quenched level 5 , and the lattice data in the chi- 
ral limit is well reproduced by 



M(p2) = Afo/{l + {p/p)-*} 



(3) 



with A/o=260 MeV, p=870MeV and 7=3.04. The 
infrared quark mass M(0) = Mq ~ 260MeV 
seems consistent with the constituent quark mass 
in the quark model. Using this lattice result 
of M(p^), the pion decay constant is calculated 



2 



as f„ ~ 87 MeV with the Pagels-Stokar for- 
mula and the quark condensate is obtained 
as {qq)A=iGcV — — (220MeV)'^. These quantities 
related to DCSB seem consistent with the stan- 
dard values. 



to be approximated as Z^p'^) — 1, which reduces 
the SD equation to 

M( -xr f ^^(g') Kjip-q)^) 



3. The Schwinger-Dyson Equation 

In this section, we formulate the SD equation 
for quarks in the chiral limit in the Landau gauge. 
By taking the trace Eq.(Q, one finds 



M{p^) 



Z{q^) 



A + M{q^) 



T,{p,q)}D^,{q) (4) 



with q = p — q. For the quark-gluon vertex, we 
assume the chiral-preserving vector-type vertex. 



(5) 



which keeps the chiral symmetry properly. (In 
contrast, to be strict, the Higashijima-Miransky 
approximation jll2| explicitly breaks the chiral 
symmetry in the formalism.) Then, one obtains 



M(p2) 



Z{q^)M{q^) _ 
\ q^ + M^{q^) 



r{q')D,,iq) (6) 



with Cp =4/3 being the color factor for quarks. 
In the Landau gauge, the Euclidean gluon prop- 
agator is generally expressed by 



Df,^[p ) = — 5- I - 



(7) 



where we refer to d{p^) as the gluon polarization 
factor. Therefore, Eq.® is rewritten as 



Here, we define the kernel function 
K{p') ^ g'r{p')d{p') 



(8) 



(9) 



as the product of the quark-gluon vertex r(p^) 
and the gluon polarization factor d{p^). Then, 
the SD equation is expressed as 



Z{q^)M{q^) mp-qf) 
' q^ + NP{q^) [p-qY ' ^^^> 



In the Landau gauge, the quark wave-function 
renormalization is not so significant and seems 



4. Extraction of the Kernel Function in the 
SD Equation from Lattice QCD 

In this section, we extract the kernel function 
K{p^) EE g'^V{p^)d{p^) in the SD equation ^ 
using the quark mass function M{p^) obtained 
in lattice QCD. By shifting the integral variable 
from q to q = p — q, we rewrite Ea. pi|l as 



M{p^) = 3Cf 



M{{p-q) 



Km 



,^{p^qY + M\{p-qY) Ip 
Therefore, we obtain 



(12) 



M(p2 



3Cf 



dq'Q{pA)K{q'), 



J, 

where Q{p, q) is defined with M{p'^) as 



(13) 



Q[p,q) = [ dOsuY 
Jo 



M{p^ + q^ — 2pq cos t 



p2 ^ q2 _ 2pq cos 9 + AI^ [p^ -f — 2pq cos 6) 



(14) 



Since the quark mass function M{p'^) is given by 
Eq.lO in lattice QCD, we can calculate the kernel 
function Klp^) from Ea. fT^ . 

As shown in Fig.l, we numerically obtained the 
kernel function K{p^) = g^V{p^)d{p^) extracted 
from the lattice QCD result for the quark propa- 
gator in the Landau gauge. 

As remarkable features, we find "infrared van- 
ishing" and "intermediate enhancement" in the 
kernel function K{p^) in the SD equation. In fact, 
K{p^) seems consistent with zero in the very in- 
frared region as 



Kip" ^ 0) ^ 0, 



(15) 



while K{p^) exhibits a large enhancement in the 
intermediate-energy region around p 0.5GeV. 

These tendencies of infrared vanishing and in- 
termediate enhancement in the kernel function 
K{p'^) = g^V{p^)d{p^) are observed also in the 
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[GeV'] 

Figure 1. The kernel function in the SD equation, 
K^p^) = g^r{p'^)d{p^), extracted from the lattice 
QCD result of the quark propagator in the Lan- 
dau gauge. The calculated data are denoted by 
the square symbols, and the solid curve denotes 
a fit function for them. In the very infrared re- 
gion, the kernel seems consistent with zero. As 
remarkable features, K{p'^) exhibits infrared van- 
ishing and intermediate enhancement. 



direct lattice-QCD measurement for the polariza- 
tion factor d{p'^) in the gluon propagator in the 
Landau gauge j7 . 

Note here that the usage of the perturbative 
gluon propagator and the one-loop running cou- 
pling grun(p^) corresponds to if(p^) = 3?un(P^)> 
which largely differs from the present result based 
on lattice QCD both in the infrared and in the 
intermediate-energy regions. In fact, the simple 
version of the SD equation using the perturba- 
tive gluon propagator and the one-loop running 
coupling would be too crude for the quantitative 
study of QCD. 



5. Chiral Symmetry at Finite Temperature 

Finally, we demonstrate a simple application of 
the LQCD-based SD equation to chiral symmetry 
restoration in finite-temperature QCD. 

At a finite temperature T, the field variables 
obey the (anti-)periodic boundary condition in 
the imaginary-time direction, which leads to the 
SD equation for the thermal quark mass A/„(p^) 



of the Matsubara frequency w„ 



{2n 



DttT as 

■'nm ' ^ 



m.q 



M2 (q2) t^2 



q2 



,(16) 



with f = 
and q = p — q. 



Using the kernel function K{p^) obtained in the 
previous section, we solve Eq. for the thermal 
quark mass M„(p^). Figure 2 shows the prelim- 
inary result for the thermal infrared quark mass 
Afo(p^ = 0) plotted against the temperature T. 
We thus find chiral symmetry restoration at a 
critical temperature Tc ~ lOOMeV. 
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Figure 2. The thermal infrared quark mass 
Mo(p^ = 0) plotted against the temperature T. 
The critical temperature is found to be about 
lOOMeV. 
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